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Abstract 

We study thermal radiation outside equilibrium. The situation we consider consists of two bodies 
emitting photons at two different temperatures. We show that the system evolves to a stationary 
state characterized by an energy current which satisfies a law similiar to the Stefan-Boltzmann law. 
The magnitude of this current depends on the temperature of the emitters expressed through the 
difference of the fourth power of these tempereatures.The results obtained show how the classical 
laws governing the thermal radiation at equlibrium can be generalized away from equilibrium 
situations. 
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I. INTRODUCTION 



Thermal radiation at equilibrium was studied by Planck [1| by using equilibrium ther- 
modynamic concepts. The thermal properties of the gas of photons are well-known. One 
of them, the Stefan-Boltzmann law ^ gives the value of the energy flux in terms of the 
temperature of the emitter through a power law crT^. 

However, in many instances frequently found in nano-systems, the radiation is not in 
equilibrium due to the presence of thermal sources or temperature gradients. This is what 
happens in nanostructures such as solar cells and thermophotovoltaic devices la, 41, even in 
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cancer therapies |5|, 16|], just to cite some exemples. In these situations, the classical scheme 
is no longer applicable and it then becomes necessary to employ a nonequilibrium theory. A 
first attempt to describe non-equilibrium radiation could be performed via nonequilibrium 
thermodynamics Nevertheless, some of the laws governing the behavior of thermal 
radiation are non-linear laws whose derivation is beyond the scope of this theory which 
provides only linear relationships between fluxes and forces. 

In this article, we will show that this limitation can be overcome if mesoscopic nonequi- 
librium thermodynamics jo], Ql is used. This can be done by performing a description in 



terms of an internal variable 111 



12| -the momenta of the photons- and by assuming local 



equilibrium in phase-space. In this way it is possible to analyze the underlying activated 
process for which photons are emitted. The photon energy current does obtained leads to 
the nonequilibrium Stefan-Boltzmann law. 

This article is organized as follows. In Section 2, we study the photon gas outside equi- 
librium. Starting from the Gibbs entropy postulate, we derive the entropy production and 
from this the current of photons. In Section 3, we analyze the steady state obtained when 
photons are emitted by two bodies at two different temperatures. By using mesoscopic 
nonequilibrium thermodynamics we derive the equivalent to the Stefan-Boltzmann law in 
this nonequilibrium situation. In the conclusion, we discuss some perspectives of the results 
obtained. 



2 



II. THE NONEQUILIBRIUM PHOTON GAS 



Let us consider a gas of photons distributed through the law n{T,t) which is the prob- 
abihty density defined in the single-particle phase-space T = (p,x), where p, x are the 
momentum and position of a photon, respectively. According to the principle of the conser- 
vation of probabihty, we assume that n{T,t) satisfy the continuity equation 

|^(r,t) = -|t • J (r,t) (1) 

where d/dT = (d/dpjV), V = d/dx. The continuity equation ([1]) defines the probability 
current J (r,t) = (J^,, Jp) which must be determined by means of the methods of nonequi- 
librium thermodynamics. 

According to the Gibbs entropy postulate, our thermodynamic analysis of the nonequi- 
librium gas is based in the assumption of the density functional 



Sit) 



-ks [ n{T,t)ln^^^dT+So (2) 

J fT'eqi^ ) 



as the nonequilibrium entropy of the system |8l-ll0| . Here, So is the equilibrium entropy of 
the gas plus the thermal bath and neq{T) is the equilibrium probability density function. 
By taking variations of Eq. ([2]) we obtain 

5S = -ks [ Sn{T,t) In ^^^dT , (3) 
where once one introduces the nonequilibrium chemical potential 

Mr,t) = A;^rin^, (4) 



Eq. ([3]) can be written 



5S = - I ^5n{T,t)dT . (5) 



Here, Eq. ([5]) which is the Gibb's equation of thermodynamics formulated in the phase- 
space, illustrates the physical meaning of the nonequilibrium chemical potential (jl]). Since 
—T6S = 6F, with F the nonequihbrium free energy, Eq. ([5]) leads to 

SF = J fi{T,t)6n(T,t)dr . (6) 
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On the other hand, the equihbrium chemical potential of photons is zero, thus no reference 
value /io is needed in Eq. (j4]). From Eqs. ([1]) and ([5]) we obtain the entropy production 

f^-/.(r.,.A^.r, 

as the product of a thermodynamics current and the conjugated thermodynamic force 
d/dT{fi{T,t)/T) [8]. As usual in nonequilibrium thermodynamic js], from Eq. ([7]) we derive 
the phenomenological law 

j(r,) = -L(r).A^i(M . (8) 

where L(r) is the matrix of phenomenological coefficients. Hence, Eq. ([8]) enables us to 
write the entropy production Eq. ([7]) as a bilinear form 

3S_ fa ,(r,t) j^i^j 9 MM,r . (9) 



dt J dT T ^ ' dT T 

In terms of the diffusion matrix D(r) = /i;_bM, where M(r) = L/n is the mobility, Eq. ([8]) 
can be written in a more convenient form 

j(r,t) = -D(r)- A^r,t), (lo) 

which constitutes the Pick's law of diffusion formulated in the single-particle phase- 
space. 



III. STATIONARY STATE AND THE NONEQUILIBRIUM STEFAN- 
BOLTZMANN LAW 

In this section we will study the heat exchange by thermal radiation between two bodies at 
different temperatures. This process comes into play in nanostructures such as solar cells and 
thermophotovoltaic devices, just to cite some exemples. Applications of thermophotovoltaic 
devices range from hybrid electric vehicles to power sources for microelectronic systems j^. 
Therefore, as mentioned, the recent developments of nanotechnology made the study of heat 
exchange by thermal radiation an object of growing interest jisl . 

Hence, to undertake this study, we assume that the dynamics of the photons is the result 
of two simultaneous processes: emission and absortion of cold photons at Tc and emission 
and absorption of hot photons at Th, this is illustrated by the Figure (see Fig. (1)) 



4 




FIG. 1: (Color online) Schematic illustration of the radiation exchanged between two materials 
maintained at different temperatures, Th and Tc, separated by a distance L. 

Since the photons do not interact among themselves we assume that the system is homo- 
geneous, r — )■ p, and the diffusion matrix reduces to a scalar -D(r), the diffusion coefficient. 
Additionally, if there are only hot and cold photons 



J (p,t) = Jc(t)5(p - Pc) + J/f(t)5(p - Vh) 



(11) 



i.e. the system reaches a state of quasi equilibrium. Thus, integration of Eq. flTU]) taking 
into account Eq. ffTTl) leads to 



Jc{t) , Juit) 



[n(pc,t) - n(p//,t)] , 



(12) 



Dc Dh 

with J{t) = u • J(t) and u is the unit vector normal to the walls. Whence, introducing the 
net current J{t) defined through 

J{t) Jcit) , Juit) 



aDc'Dn D. 



c 



D 



H 



or equivalently 



J{t) = aDnJcit) + aDcJH{t) 



(13) 



(14) 
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we can rewrite Eq. (fT2|) 



J{t) = -aDcDu Hvc.t) - n{^H,t)] , (15) 

where a is an effective parameter accounting for the dimensionahty. Here, term by term 
comparison of Eqs. ( IT^ and (fT5|) leads to the identification 

Jc{t) = -aDcn{pc,t), (16) 

Jnit) = aDHn{pH,t) . (17) 

Therefore, 

- DnJcit) = DnDcnipct) (18) 

represents the fraction of photons absorved at the hot wall from the fraction Jc{t) of photons 
emited at the cold wall. In the same way, 

DcMt) = DcDHnipH,t) (19) 

represents the fraction of photons absorved at the cold wall from the fraction Jnit) of 
photons emited at the hot wall. 

To better illustrate the physics under the nonequilibrium process involving the photons, we 
introduce the affinity 

A = fic- /^iT, (20) 

where fic = fJ'ipct) and fin = fJ'{PH,t), which with the help of Eq. (jl]), enables us to rewrite 
Eq. ([I5D as [ill 

This result shows the nonlinearity inherent to the process involving the radiation since this 
corresponds to a heat current which does not satisfies the Fourier law. 

In general, the diffusion coefficient might depend on the frequency. However, one can 
introduce a cut-off limit. Ay = cH/ksT, the thermal wavelength of a photon, which marks 
this dependence. So, for length scales L ^ Xt {i.e. high frequencies) we can treat the 
photons as point particles, hence in this case we deal with the free diffusion of point particles 
for which the diffusion coefficient is constant. On the other hand, when L < Xj- [i.e. low 
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J(t) = aDcDu exp 



1 — exp 



A 



exp 



frequencies) we deal with the problem of casre-diffusion 16 1, therefore the diffusion coefficient 

nn 

must depend on the ratio \t/L or equivalently on the frequency |17l. Il8l|. 

In the stationary state, the current of photons per unit of volume in phase-space is given 

by 

Jst{uj) = aDc{oo)DH{oo) [n{uj, Th) - n{u, Tc)] , (22) 

where 

nKr)=2^, (23) 



with h being the Planck constant and N{ijj^ T) the averaged number of particles in a ele- 
mentary cell of the phase-space given by the Planck distribution {2] 



^<"-^> = exp(ft.Ar)-r (2^' 

In Eq. fl23|) the factor 2 comes from the polarization of the photons. Here, the diffusion 
constant plays the role of the effective cross section. 

By multiplying Eq. fl22|) by hw -the energy of a photon- and integrating in momentum 
space we obtain the heat current 

Q = J fiojJsti(^)dp, (25) 

where p = {}iijj/c)flp. In the case L Xt , and in an ideal situation, we can take 
Dc = Dh = D = 1. Thus, with Eq. dMD, Eq. ^ leads to 

Q = ^ J dudnpu'Muj) = (y{Ttj- T^) , (26) 

where we have taken a = c/4, with c the velocity of the radiation and a = n^k^/GOH^c^ 
being the Stefan constantn [2]. At equilibrium Th = Tc, therefore Q = 0. 

For a s-dimensional system, the phase-space is made up of s coordinates and s conjugated 
momenta. Hence, the volume of an elementary cell of this phase-space is and now 

ni.,T) = 2^. (27) 

In this case, by integrating over a hypersphere in momentum-space and assuming that now 
a is a function of s, a{s), it this possible to obtain 

Q = a{s) (T^+i - T^+i) (28) 

which generalizes results previously derived in Ref. 21 1 



IV. CONCLUSIONS 



In this paper we have performed a thermodynamic analysis of the radiative heat ex- 
change between two bodies at different temperatures separated by a certain distance L. 
We have shown how the Stefan-Bolztmann law can be generalized to nonequilibrium situa- 
tions as those in which thermal radiation is composed by photons emitted at two different 
temperatures and thus having different momenta. In the framework of mesoscopic nonequi- 
librium thermodynamics based on the Gibb's entropy postulate, we have obtained the Gibb's 
equation of thermodynamics which here, describes the local equilibrium in the phase-space. 
Then, by means of the usual procedure of nonequilibrium thermodynamics, from the Gibb's 
equation we have derived the Pick's law for the diffusion of photons valid at short time 
scales. 

The exchange of energy is assumed to be due to activated processes related to reaction rate 
currents which are derivable in the framework of our thermodynamic theory. These currents 
provide us with the rate of absortion and emission of photons. Since in the stationary state 
both hot and cold photons are in local equilibrium, their rates are proportional to the Planck 
distribution. Pinally, the net current of heat is given through a balance of these rate currents 
after integration over frequencies, constituting the nonequilibrium Stefan-Boltzmann law. 

Systems outside equilibrium exhibit peculiar features not observed in an equilibrium 
state [7]. The results obtained show how non-equilibirum phenomena taking place in gases 
composed of quasi-particles, governed by non-linear laws, can be analyzed by means of the 
methods of nonequilibrium thermodynamics applied to the mesoscale. These methods have 



been shown to be very useful in the study of activated processes of different natures 
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